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Abstract 

We study the speed of propagation of fronts for the scalar reaction-diffusion 

equation ut = Uxx + f{u) with /(O) = /(I) = 0. We give a new integral 

variational principle for the speed of the fronts joining the state u = 1 to 

u = 0. No assumptions are made on the reaction term f{u) other than those 

needed to guarantee the existence of the front. Therefore our results apply to 

the classical case /> in (0, 1), to the bistable case and to cases in which / 

has more than one internal zero in (0, 1). 
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Typeset using REVT^ 



The one dimensional reaction diffusion equation 

ut = u,. + fiu) with /(0) = /(l)=0 (1) 

with f{u) G C^[0,1] has been the subject of much study as it models diverse phenomena 
in biology, population dynamics, chemical physics, combustion and others Not only 

is it itself of interest, but based on the rigourous results available for this equation, diverse 
methods applicable to pattern forming systems have been developed [|]-|^. In applications, 
the reaction term f{u) obeys additional requirements depending on the phenomenon being 
modelled. Three types of nonlinearities appear to be generic, and are shown in Fig. 1. 



Type A, for which / > in (0, 1) is the class to which the classical case of Fisher and 
Kolmogorov, Petrovskii and Piskounov (KPP) |ll| belongs. Type B, usually referred to as 
the combustion case satisfies / = on (0, a) and / > on (a, 1), while finally, type C, called 
the bistable case, satisfies f{u) < for m in (0, a), / > on (a, 1) with /o f{u) du > 0. More 
general cases for /, namely cases in which / has more than one internal zero have also been 

studied mm. 



The time evolution of an initial condition u{x, 0) has been studied for all the cases 
mentioned. It was proved that for suitable initial conditions the disturbance evolves 



into a monotonic travelling front u = q{x — ct) joining the stable state m = 1 to m = 0. 
In case A there is continuum of values of c for which a monotonic front exists, the system 
evolves into the front of minimal speed. In cases B and C there is a single isolated value 
of the speed for which the front exists. In these two last cases there are threshold effects, 
necessary conditions for the evolution of the system into the front have been established 



as well [lT2| , p^ . The same is true for reaction terms with more than one internal zero [[12 
The problem which interests us here is the determination of the asymptotic speed of the 
front. There have been numerous studies of this problem, a very complete review is given 
in H]. For reaction terms of type A which in addition satisfy /'(O) > f{u)/u the speed is 



given by ||n| c = ckpp = 2y/'(0). Another reaction term of type A, that of a function 
/ approaching a Dirac function at u = 1 was studied related to combustion by Zeldovich 



and Frank-Kamenetskii. They showed [0 that the speed tends io czfk = y'^ f{u) du. 
For an arbitrary reaction function f{u) a local variational principle of the minimax type 
exists 0Jl^ . For reaction terms of type A and B we have shown that an integral variational 
principle of the Rayleigh Ritz type exists |[17|JT8[| . Recently an interesting conjecture [jl9l has 



been put forward for a restricted class of reaction functions. 

The purpose of this article to show that the speed of the front joining the state u = 1 to 
M = derives from an integral variational principle without any restriction on / other than 
those needed to guarantee the existence of the front. The derivation follows an approach 
similar to the one used to obtain the principle valid only for for positive reaction terms. 
However, this new principle which is valid for all cases, is not related, nor equivalent, to the 
previous one. 

It is known ||l4| that for a function of type A, B, or C there exists a strictly decreasing 
front u = q{x — ct) joining m = 1 to m = for some c > 0. The front satisfies qzz + cqz + f{q) = 
0, hmM^^_oo = 1, limw^-^oo = 0, where z = x — ct. Following the usual procedure, since the 
front is monotonic, we define p{q) = —dq/dz, where the minus sign is included so that p is 
positive. One finds that the monotonic fronts are solutions of 

P{q) ^ - cpiq) + f{q) = 0, (2a) 

with 

p{0) = 0, p(l)=0, p>0 in (0,1). (2b) 

The derivation follows in a simple way from equation (^). Let g be any positive function in 
(0,1) such that h = —dg/dq > 0. Multiplying equation (|2a| ) by g{q) and integrating between 
g = and g = 1 we obtain, after integration by parts, 

rl rl rl I 

/ fgdq = c pgdq- -hp^ dq. (3) 
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However, since p, g and h are positive, for fixed q, the function 

0(p) = cpg - hip^ 



has a maximum at 



Pma. - (4) 



SO 

(j){p) < c' 
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at each value of g. It follows then that 

c^>2 4^- (5) 

This bound on the speed is valid for any / for which a monotonic front exists. To show 
that this is a variational principle we must show that there exists a function g at which the 
equality holds. From Eq. (|D we see that 

9 

c— = p 
h 

which can be integrated. The maximizing g is given by 

(6) 



g = expl- f -dq] 



with < go < 1- Evidently g is positive, monotonic decreasing and moreover ^(1) = 0. Near 
q = g diverges. We must ensure that the integrals in Eq. (^) exist. To verify this we recall 



14 1 that in the three cases, A, B and C, the front approaches g = exponentially, therefore. 



Then we obtain 

9{(l) 



qc/m. 



near zero and fg and g'^/h diverge at most as g^ (c/m)_ rjj^g integrals in Eq.(|^) exist if 
m/c > 1/2. This condition is always satisfied when /'(O) < 0, that is in cases B and C. 



In case A this condition is satisfied provided that c > 2y/'(0), which is the KPP case. 
However, choosing as a trial function g{q) = a{2 — a)-u"~^ with < a < 1 one can check 
that in the limit a —>■ 0, —>■ 4/'(0). Therefore, including all cases our main result is 



where the supremum is taken over all positive decreasing functions g in (0, 1) for which the 
integrals exist. Moreover there is always a maximizing g in cases B and C, while for case A 
there is a maximizing g whenever c> ckpp- 

While we considered the case of decreasing fronts, similar results can be derived for 
increasing fronts and for the case of density dependent diffusion following the same approach. 

As an example we may apply the above result to the Nagumo equation which corresponds 
to a reaction term of the form 

f(u) = u{l-u){u-a) with < a < 1/2 

This reaction term is of the bistable type. For this equation the solution to Eq.(^) is known, 
it is given by p{q) = — q) and the speed is given by 

c = —pz — a\pi. 
V2 

To exhibit in this solvable case that the exact speed can be obtained from the variational 
principle choose as a trial function 



\ l-2a 

1 -g\ 



The integrals can be performed easily. We obtain 



and 



^ ^ > ^ (1 -2a)r(4) 



1 , , (l-2a)r(l + 2a)r(3-2a) 



/o ' ^ ^ 4r(4) 
so that = (1 — laf' j1 which is the exact value. For other non solvable cases, it is a simple 

matter to obtain accurate values for the speed using standard variational techniques. 
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